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Abstract 

As a part of a thorough work of excitation functions of deuteron induced reactions, experimental cross-sections of 
i85,i83m,i83g,i82Qg I88,i86,i84m,i84|;,i83j^g activation products on ""'Re Were measured up to 40 MeV for the first time 
with the activation method using a stacked foil irradiation technique and high resolution -spectrometry. Comparison 
with the former results of other laboratories and with the predictions of the ALICE-IPPE and EMPIRE-3 model 
codes, modified for improved calculations for deuteron reactions, and with data in the TENDL-2011 library are also 
presented. Thick target yields were given deduced from our experimental cross-sections and compared with the 
few literature values. For practical applications (thin layer activation) also activity versus depth distributions were 
calculated for selected isotopes. 
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1. Introduction 

Integral excitation functions for the production of 
residual nuclides through light charged particle acti- 
vation constitutes basic data for different applications. 
To meet the requirements of different practical appli- 
cations we started to establish an experimental activa- 
tion database some years ago by performing new ex- 
periments and a systematical survey of existing data of 
deuteron induced cross-sections up to 50 MeV We 
present here new results on deuteron activation cross- 
section data on rhenium. No earlier experimental cross- 
section data were found in the literature. Only a few ex- 
perimental thick target yields are available at 22 MeV 
reported by Dmitriev et al. [21. Rhenium (atomic 
mass: 186.207 g/mol, melting point 3180.0 "C, density 
2 1 .02 gjcnr' ) is a very heavy (dense), high melting point 
metal. Rhenium is resistant to heat, wear and coiTosion. 
It is able to maintain ductility between absolute zero and 
its melting point 3180 °C. This makes it valuable as 
an alloying agent for hardening metal components, as 
well as for use under extreme temperatures. Rhenium 
is used in many fields (catalyst, electric components. 
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coatings, etc.). Presently it is mostly used as an addi- 
tive in super-alloys for aviation technology. Activation 
cross-sections of proton and deuteron induced reactions 
are important therefore for estimation of activation dose 
and for use in thin layer activation technology. A new 
application is the production of medical radioisotopes. 
Some alternate production routes were investigated in 
detail earher li3n6J. 

2. Experimental method and data evaluation 

The experimental method was practically the same 
as that we applied before for charged particle in- 
duced nuclear reactions measurements |7|. Activation 
method using stacked foils was applied for measur- 
ing the whole excitation functions. Two foil-stacks 
were irradiated at 21 MeV (VUB) and one stack us- 
ing 40 MeV (CYRIC) deuteron energy. In our experi- 
ments high-purity Mo(52.5)Re(47.5) alloy foils (Good- 
fellow > 99.98%, thickness 50 /im) were stacked to- 
gether with Ti monitor foils (10.9 m) at VUB and Al 
degrader/monitor foils (100 m) at CYRIC. The "long" 
target stack at CYRIC also contained Rh (12.3 fim). Ho 
(25 urn) Au (107 iim) Yb (25 ftm), Ni(2)Mn(12)Cu(86) 
alloy (25 ^m) foils for simultaneous measurement of ac- 
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tivation cross-sections on these elements. The Ti and 
Al monitors foils were used to catch recoiled radioac- 
tive products as well as for monitoring the beam inten- 
sity and energy by comparing the excitation function of 
the "'"Ti{d, x)48y reaction at VUB and ^'^ Al{d, xf^'^'^Na 
reactions at CYRIC over the entire energy range. To 
measure the nominal beam current the stacks were ir- 
radiated in a target holder, which served as a Faraday 
cup and was equipped with collimator (beam diameter 
on target is 5 mm) and a secondary electron suppressor 
unit. Irradiations were performed using constant beam 
current of 160 nA for 60 min (VUB) and 24 nA for 30 
min (CYRIC) respectively. The number of the produced 
radionuclides was determined through their decay. For 
measurement of the activities in the target and the mon- 
itor foils HPGe -ray spectrometers were used. Mea- 
surements of the induced activity started in Brussels a 
few hours after BOB (End of Bombardment), while in 
Sendai a cooling time of at least a day was set in order to 
get rid the high activity short-lived isotopes. Chemical 
separation was not necessary and the measurements in 
Brussels were performed repeatedly several times (up to 
months) after BOB, ensuring accurate determination of 
both short-lived and long-lived activation products (the 
possibility to measure short-lived activities has been lost 
in the Sendai experiment). The samples were mea- 
sured at appropriate detector- sample distances (70 cm 
to 5 cm) to ensure low dead-times and to decrease pile- 
ups. The decay- and other nuclear data were collected 
from the NUDAT2 data-base [8| and are given in Table 
1 . The Q- values refer to formation of the ground states 
and taken from fW\. The cross-sections were calculated 
from the well-known activation formula with measured 
activity, particle flux and number of target nuclei. A 
number of radionuclides come as a result of cumulative 
processes, as decay of parent nuclides gives contribu- 
tion to the production. The particular situation for the 
individual nuclides will be explained case by case. El- 
emental rhenium has two stable isotopes with the fol- 
lowing abundances: ^^^Re(37A0%) and Re(62.60%). 
The absence of ^^^Re (not stable) in principle allows in 
many investigated cases and in special energy ranges to 
obtain pure isotopic cross-sections, required for com- 
parison with model calculations. As in practical appli- 
cations mostly natural rhenium is used, we have de- 
duced so-called elemental cross-sections, considering 
the target as monoisotopic. The number of the target 
nuclei was determined from the surface density of the 
foils obtained by precise measurement of the weight 
and the surface of the metal sheets used before cutting 
into pieces. The thickness of each individual target and 
monitor foil was checked for control after cutting. Ac- 



cording to the previous practice the beam current and 
the incident beam energy were determined by using 
the excitation function of the re-measured monitor re- 
actions IfTOl [TTI . The energy loss versus depth function 
in the stack was calculated 1 12 1 and refined using the re- 
measured excitation functions of the monitor reactions 
as given in lfT3l . The latest available decay data were 
used for our calculations |8 1. The cross-section curve of 
the simultaneously measured monitor reaction is com- 
pared with the recommended data as described in our 
previous works [14|. Acceptably good agreement was 
found after moderate adjustment on the beam current 
(10% compared to the Faraday cup measurement) and 
on the incident beam energy (0.5 MeV). The resulting 
uncertainty of the cross-sections contains the individ- 
ual uncertainties of the processes contributing linearly 
to the final result accordingly to the well-accepted sum- 
mation rules [15J. The uncertainties of the non-linear 
processes like half-life, irradiation time, measuring time 
were neglected in the reported uncertainties of the re- 
sulting cross-sections. The uncertainty of the energy in 
the targets was estimated from the uncertainty on the en- 
ergy of the primary beam, the uncertainty in the thick- 
ness and uniformity of all foils, the beam straggling and 
the well-known cumulative eff'ects. 

3. Tehoretical calculations 

In order to avoid large mistakes in evaluations of 
the measurements, to analyze and understand contribu- 
tions of individual reactions and to check the predic- 
tivity of recently used model codes, the ALICE-IPPE 
L16J, EMPIRE nil and TALYS LIU codes are used in 
our works. While ALICE-IPPE calculates only the to- 
tal cross-section of the concerned reaction channels, the 
TALYS and EMPIRE codes permit to calculate a popu- 
lation of different low-lying levels and can thus estimate 
the isomeric ratios for these levels. Independent data for 
isomers with ALICE-D code was obtained by using the 
isomeric ratios calculated with EMPIRE. During the re- 
cent analyses of several deuteron induced reactions we 
were confronted with poor description of the measured 
cross-sections due to the high complexity of the interac- 
tions. The interaction of deuterons with the target nuclei 
proceeds largely through direct reaction (DR) processes 
for deuteron energies below and around the Coulomb 
barrier, while with increasing incident energy, reaction 
mechanisms such as pre-equilibrium emission (PE) or 
evaporation from the fully equilibrated compound nu- 
cleus (CN) also become important. The breakup mech- 
anism is responsible for the enhancement of a large vari- 
ety of reactions along the whole incident-energy range. 
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and thus its contribution to the activation cross-sections 
has to be expHcitly taken into account |19|. There 
are different possibilities to get better agreement. To 
get better descriptions in our approach modified codes 
named ALICE-D and EMPIRE-D were developed at 
IPPE [14. .20.1 . in which a simulation of direct (d,p) and 
(d,t) transitions with general relations for nucleon trans- 
fer probability in the continuum |21 1 is included using 
an energy dependent enhancement factor for the corre- 
sponding transitions Il22l . These updated codes were 
used to analyze the present experimental results. The 
parameters for the optical model, level densities and 
pre-equilibrium contributions were taken as described 
in f23l. In the figures we also present the theoretical re- 
sults presented in TENDL 2011 data library L24il . The 
TALYS and TENDL calculations do include a breakup 
component in all (d,n) and (d,p) reaction channels, ac- 
cording to our earlier experiences, its strength does not 
show enough enhancement. 



4. Results 

The new experimental cross-sections are shown in 
Figs. 1-9 together with results of the theoretical calcu- 
lations. We present the results of different irradiations 
separately to show the possible systematic deviations 
and errors. The numerical data essential in further data 
evaluation are collected in Table 2. The main contribut- 
ing processes are shown in Table 1, together with the 
reaction Q- values. 

4.1. Cross-sections 

4.1.1. Production of osmium isotopes 

The measured excitation functions came in all cases 
from a combination of (d,xn) reactions on the "*^/?e and 
'^''Re. 

""'Re{d,xy^^Os 

The main contributing processes are the direct pro- 
duction via the '^^^Re{d,2n) and Re{d,Ari) reactions. 
The experimental and theoretical excitation functions 
are shown in Fig. 1 . The results of different irradiations 
show acceptably good agreement in the overlapping en- 
ergy region. The theoretical predictions differ by factor 
of two. 

""'Re(d, and ""'Reid, x^^^Ws 

The Q-value of the two contributing reactions 
'^^/?e(J,4n) and ^^^Re(d,6n) are collected in Table 1. 
We present independent cross-sections for production 
of the metastable state (Fig. 2) and for the ground 
state (Fig. 3) obtained after a small correction (low cr„. 



"'"Re(d,x}"»Os 
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Figure 1: Experimental and theoretical excitation functions for 
"'"Re{d,x)'^^Os 



IT=15 %) from the decay of metastable state. Theo- 
ries in average reproduce the shape and the magnitude 
of the experimental data, but in the investigated energy 
range there are large differences in the magnitude of the 
different model predictions. 




Figure 2: Experimental and theoretical excitation functions for 

""'Reid, xy^^Os 

Only in a few foils in the beginning of the stack we 
can detect the y-lines from the decay of the 1820i. The 
few obtained experimental data points fit well to the the- 
oretical excitation functions (see Fig. 4). In the inves- 
tigated energy range only the ^^^Re(d,5n) contribute to 
the production of ^^^Os. 
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Table 1 : Decay characteristics of tlie investigated activation products 



Nuclide 


Half-life 


E (keV) 




Contributing reactions 


Q-vatue (keV) 












^019.71 


e: 100% 


93.6 d 


646.116 


78 


^^"^ Re(d,4n) 


-17555,9 


— la'^niQ 

IT; 15% 






1034.86 
1101.92 


6,02 Tf,. 


-19479,6 


e: 85% 


170.715 keV 


9.9 h 


1107.9 


22,4 ^^''Reid,(,n) 


-33015,8 




■ 


114.47 


20.6 


185 D^t^ A„\ 


-19308.9 


e: 100% 


381.763 


89.6 




-32845.1 


182/1 


21.84h 


130.80 


3.30 


"*^Se(ii,5>i) 


-26433,9 


e: 100% 


180.20 


34.1 


l*'se(ii,7n) 


-39970.1 


188 D 












P~ : 100% 


17.004 h 


155.041 


15.61 


1 87 n , > » 


3647.184 










l*5«e(<(.p) 




/T: 92.53% 








3954.794 


e: 7.47% 


3,7183 d 


137.157 


9.47 




-9581.4 






104 7395 








IT'74 5% 




2K 547 


6^56 
95 


«,!(</, p2n) 


-10079 4 






"'5^845 


10 8 


^^'^ Re(d p4n) 


— -23545.6 







1 11.2174 


17.2 










792.067 
894.760 


37.7 
15.7 


^^^Re(d. pin) 


-9891.4 


e: 100% 


35.4 d 


903.282 


38.1 


^^'^Re(d.p4ti) 


-23427.59 






162.3266 


23,3 


^^^Re(d,p3n) 


-16378.32 


e: 100% 


70.0 d 


291.7282 


3.05 


^^'^Re(_d,p5n) 


-29914.52 






100.12 


14,4 










152.43 


7,0 










470.26* 


2,02 










894.85* 


2,11 






182m/jg 
e: 100% 




1121.4 
1189.2 


32,0 
15,1 


^^^Re(d.p4n) 


-24813,0-F 


0+X keV 


12.7 h 


1221.5 


25,0 


^^'^Re(d,p6n) 


-38349,0-1- 






100.10 


16,5 










130.81* 


7,5 










169.15* 


11.4 










191.39* 


6.7 










229.32 


25.8 










286.56* 


7.1 










351.07* 


10.3 










1076.2* 


10,6 










1121.3 


22,1 










1221.4 
1231.0 


17,5 
14,9 


^^^Reid,p^) 


-24813,0 


e: 100% 


64.0 h 


1427.3* 


9.8 


187fle(d,p6n) 


-38349.0 



+ The energy of the excited state is unknown 

When complex particles are emitted instead of individual protons and neutrons the Q- values have to be decreased by the respective binding 
energies of the compound particles: np-d, +2.2 MeV; 2np-t, +8.48 MeV; n2p-3He, +7.72 MeV; 2n2p- , +28.30 MeV 
The independent y-Unes are marked with * 



4.7.2. Production of rhenium isotopes 

The measured excitation functions came in all cases 
from a combination of (d,xn) reactions on the ^^^Re and 
'**^7?e. The level schemes for ^^^Re and ^^'^Re are es- 
timated with large uncertainties for some important y- 
transitions. So, an accuracy of calculations for the cor- 
responding isomer yields cannot be very high. 

'""Re{d,xy^^Re 

The ^^^Re is produced only through ^^^Re{d,p) re- 
action. The ^^^Re could also be produced from 
^^''Re through (n,y ) reaction by neutrons coming from 
deuteron break-up. This effect was excluded, because 
out of the target foils referred in the experimental 
part one more foil of the same material was inserted 
at the end of the stack beyond the threshold of the 
""'Re{d,xy'^^Re reaction and it did not contain ^^^Re. 
It means that the neutron induced production is negU- 



gible. The measured and calculated excitation functions 
are shown in Fig. 5. There is a strong underestimation 
of the measured cross-sections in case of TENDL 201 1 . 
The agreement with the maximum value of the updated 
version ALICE-D and EMPIR-D is better, but stiU with 
underestimation. In case of EMPIRE-D the reproduc- 
tion of the shape of the high energy tail of experimental 
excitation function is also poor. 
"'"Re{d,xy^^«Re 

The ^^^^Re is produced through the ^^^Re{d,p) and 
^^^Re(d, p2n) reactions. As discussed in in the previous 
section the neutron induced production is also negligi- 
ble here because of the lack of ^^^^Re in the last foil. 
There is a strong underestimation of the measured cross- 
sections in case of all models for (d,p) reaction (Fig. 
6). The situation is the same in the case of (d,p2n) for 
ALICE-D, EMPIRE-D and TALYS (TENDL-2011). 

'^'Re{d, xf^"'Re and xf^«Re 
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Figure 3: Experimental and theoretical excitation functions for 



Figure 4: Experimental and theoretical excitation functions for 



The experimental and theoretical cross-sections are 
shown in Figs. 7-8. The models describe well the shape, 
but the estimation of the magnitude is not so successful, 
probably due to the unclear level scheme. 

""'Re(d,xy^^Re (cum) 

The ^^^Re {T1/2 = 70. OJ) is produced directly via 
(d,pxn) reactions C^^Re(d, p3n) +'**^ Re(d,p5n) and 
by decay of the ^^^""Os (e:85%, Tj/j = 9.9h), and 
■**^«Oi(e:100%, ri/2 = 13.0/!). Our cross-sections were 
deduced after complete decay of the parent isomers, in 
such a way they represent cumulative cross-sections. 
The measured and calculated cross-sections are shown 
in Fig. 9. The ALICE-D, EMPIRE-D and TENDL- 
2011 describe approximately the shape of the experi- 
mental excitation function. In magnitude the best agree- 
ment is with ALICE-D. 

"'"Re(d,xy^^'"Re and natReid, xY^^i^Re 
The "*2m/;e = i2Jh) is produced both directly 
and by the decay of the simultaneously produced '^^Os. 
The '^^Os (7^1/2 = 2lMh) is decaying only to ^^^"'Re 
(Ti/2 = 12.7/z). From other side the ^^^«Re (64 h) 
is only produced by direct production route (no decay 
contribution from ^^^Os and no internal transition from 
the isomer. According to the theoretical estimations 
the cross-sections for the direct production of both iso- 
meric states are low in the investigated energy range 
and the Q value of ^^^Re{d, p4n) reaction is high (- 
24.8 MeV). As it was mentioned, the ^^^"'Re is pro- 
duced also through the decay of the significantly higher 
cross-section ^^^Os. The "*^'"7?e has strong common 
y-lines with the ^^^^Re (see Table 1.) and a few in- 



dependent weak y-Hnes (470.2 keV and 894.85 keV). 
We could detect these independent y-lines in our spec- 
tra only with poor statistics. We have recognized that 
the separation of the high yield contribution of the par- 
ent could be done only with very large uncertainties 
and so with low reliability. Therefore we have tried 
to find direct production cross-section of the 182gRe. 
The ^^^^Re has a few moderately strong independent y- 
lines: 130.81 keV(7.5%), 169.15 keV(l 1.4%), 191.39 
keV(6.7%), 286.56 keV(7. 1 %), 35 1 .07 ke V( 10.3%) and 
1076.2 keV(10.6%). Unfortunately due to the complex 
y-spectra caused by numerous reactions on Mo and Re 
target isotopes, the reliable identification of these y- 
lines was unsuccessful. The theoretical codes give also 
very discrepant results for these two reactions. 



5. Thick target yields and activity distribution 

In different practical applications either the total pro- 
duced activity or the distribution of that activity as a 
function of bombarding energy and/or the depth are of- 
ten more convenient to use than the microscopic cross- 
section data. The production yield is directly connected 
to the cross-section and to the composition of the used 
target. From fits to our experimental excitation func- 
tions thick target yields were calculated and are shown 
in Fig 10 in comparison with the experimental thick tar- 
get yield data in the literature. The deduced yields are 
so called physical yields, calculated for an instantaneous 
irradiation 1 25 1 . Experimental thick target yields were 
found only for 1^6) and ^^^Re ||27| in the literature 
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Table 2: Measured cross-sections for the Os and Re radioisotopes 



series 


Bombarding 


185 


0,( 




183^ 


Os 


182 







18K 


Re 


1S6 


Re 




184^ 


Re 


l^-*/?.(.um) 






±SE 


IT ±S(r 


a 1 +So- 




cr 1 ±<So- 


(T 1 +Scr 


(T 1 ±S(r 


a- 1 ±6(r 


<r 1 ±6(r 




±S<r 






mb 


mb 


mb 


mb 


mb 


mb 


mb 


mb 


mb 


1 


38.5 


0.3 


492.58 


56.49 


38.86 


4.61 


251.0 


28.3 


213.0 




24.0 


34.6 


4.4 


193.3 


22.0 


36.7 


6.5 


134.5 


15.3 


383.7 


45.3 




35.6 


0.3 


644.06 


73.50 


52.72 


6.14 


301.3 


33.9 


94.0 




10.7 


27.0 


3.9 


207.1 


23.5 


32.2 


6.0 


126.0 


14.3 


470.3 


55.2 




32.4 


0.4 


796.01 


90.15 


64.84 


7.44 


342.0 


38.5 


22.6 




2.8 


39.0 


4.8 


204.2 


23.2 


24.6 


4.7 


112.8 


12.8 


483.4 


56.0 




29.0 


0.4 


830.75 


94.08 


58.40 


6.73 


265.4 


29.9 


8.6 




1.4 


47.4 


5.7 


172.3 


19.6 


14.9 


3.7 


88.0 


10.0 


366.4 


43.2 




25.2 


0.5 


581.05 


66.11 


25.15 


3.02 


96.8 


11.0 








60.9 


7.1 


134.6 


15.4 


7.6 


3.1 


53.7 


6.1 


144.0 


18.6 




21.0 


0.5 


127.24 


16.29 
















86.1 


9.9 


97.4 


11.2 






23.3 


2.8 








16.1 


0.6 


204.27 


24.76 
















124.7 


14.2 


110.8 


12.7 






15.0 


1.8 








9.8 


0.6 


77.13 


10.22 
















96.1 


10.9 


70.3 


8.2 






3.6 


0.6 






2 


19.9 


0.3 


114.60 


13.17 
















106.4 


13.2 


100.5 


11.3 






18.7 


2.1 








18.0 


0.3 


131.14 


15.00 
















119.4 


13.6 


100.5 


11.3 






13.4 


1.6 








16.0 


0.4 


178.15 


20.67 
















136.5 


15.8 


107.4 


12.1 






10.9 


1.4 








13.8 


0.4 


244.52 


27.75 
















152.9 


17.2 


115.7 


13.0 






9.6 


1.2 








11.5 


0.4 


260.05 


29.32 
















151.3 


32.5 


122.3 


13.7 






6.4 


0.8 








8.9 


0.5 


136.23 


15.95 
















142.8 


16.1 


91.4 


10.3 






1.3 


0.5 






3 


20.2 


0.3 


61.15 


9.13 
















98.8 


17.4 


112.4 


12.6 






18.9 


2.2 








18.5 


0.4 


60.23 


1 1 .37 
















122.9 


20.3 


111.8 


12.6 






13.8 


1.7 








16.7 


0.4 


89.06 


13.57 
















144.0 


16.4 


1 16.3 


13.1 






13.1 


1.6 








14.7 


0.4 


128.14 


17.73 
















173.5 


21.0 


123.8 


13.9 






11.7 


1.5 








12.6 


0.5 


144.59 


19.10 
















183.4 


20.6 


135.7 


15.2 






10.4 


1.4 








10.1 


0.5 


100.60 


14.22 
















174.3 


19.6 


115.0 


12.9 






3.8 


0.8 









Figure 5: Experimental and theoretical excitation functions for 
'^'Re(d,xf^^Re 



Figure 6: Experimental and theoretical excitation functions for 



and are in good agreement with otir calculated results. 
Some of the investigated radionuclides (depending on 
half-Ufe, excitation functions, emitted y-hnes) produced 
in deuteron bombardment on natural rhenium are suit- 
able for thin layer activation (TLA) applications to in- 
vestigate the wear or corrosion rate of friction surfaces. 
Due to the short half-life the 186Re (3.72 d) radiotracer 
can be used for monitoring processes with quick surface 
loss, while the ^^^Os (93.6 d) is suitable for investigat- 
ing long term slow processes. The activation curves at 
optimized bombarding energies for both reactions (15 
MeV for 186Re and 13.8 MeV for ^^^Os) and 1 h and A 
irradiation and also for 10 d cooling time are shown in 
Fig. 11. 



6. Summary and conclusions 

The principal aim of this investigation was twofold: 
in the frame of a systematic study of activation cross- 
sections of the deuteron induced reactions of metals to 
measure the missing activation data on rhenium and to 
check the predictivity of the most widely used theoreti- 
cal codes. We present for the first time excitation curves 
for deuteron induced reactions on natRe targets up to 
40 MeV leading to the production of i85,i83m,i83g,i82Q^ 
jjjj^ I88,i86,i84m,i84g,i83^g radiouucUdcs. The results are 
correlating to measurements of monitor reaction over 
the whole energy range. The results of irradiations of 
difi'erent primary beam energy show good agreement in 
the overlapping energy ranges. Comparison of different 
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model codes and data libraries with experiments showed 
modest agreement, even the updated versions (modified 
for deuteron induced reactions) of the codes ALICE and 
EMPIRE are not always convincing. Although the in- 
clusion of a phenomenological enhancement factor for 
direct (d,p) reactions and explicit description of more 
compound exit channels has improved the description 
in some cases, still very large discrepancies remain in 
other cases. In average the agreement between the ex- 
perimental data and the TALYS theoretical results in 
TENDL-2011 is also not satisfactory. These new re- 
sults can have applications in accelerator technology, 
radiation protection and nuclear reaction theory devel- 
opment. The activation cross-section data of proton in- 
duced reactions in rhenium are also missing. The data 
evaluation of the proton experiments is in progress. 
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Figure 11: The and "**0.s activation curves with 15 and 13.8 

MeV bombarding energies respectively (beam current: 1 juA; irradia- 
tion time 1 h; cooling time: 10 d) 
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